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Abstract 



o ■ 

o ; 

A new hidden symmetry is exhibited in the reflection equation and related quantum integrable models. It 
q ■ is generated by a dual pair of operators {A, A*} e A subject to q~ deformed Dolan-Grady relations. Using the 
inverse scattering method, a new family of quantum integrable models is proposed. In the simplest case, the 
Q ' Hamiltonian is linear in the fundamental generators of A. For general values of q, the corresponding spectral 
I/"") • problem is quasi-exactly solvable. Several examples of two-dimensional massive/massless (boundary) integrable 
models are reconsidered in light of this approach, for which the fundamental generators of A are constructed 
explicitly and exact results are obtained. In particular, we exhibit a dynamical Askey- Wilson symmetry algebra 
, in the (boundary) sine-Gordon model and show that asymptotic (boundary) states can be expressed in terms of 
ON ' q— orthogonal polynomials. 
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o 

^ ! 1 Introduction 

The exact solution of the planar Ising model in zero magnetic field |J has provided a considerable source of 
' developments in the theory of exactly solvable systems of statistical mechanics, or quantum field theory in two 
J> . dimensions. Onsager's successful approach was originally based on the so-called Onsager algebra and its rep- 

- —i 

resentations. Using algebraic methods, Onsager derived exact results such that the largest and second largest 
eigenvalues of the transfer matrix of the model. Afterwards, the two-dimensional Ising model was reconsidered 
using the now famous free fermion (Clifford algebra) techniques [2]. The Onsager algebra itself being not neces- 
sary in this latter approach, it probably explains why it assumed only a position of an interesting curiosity in the 
following years. Despite of this, in the 1980s the Onsager algebra appeared 0E1 to be closely related with the 
quantum integrable structure discovered by Dolan and Grady in , associated with a class of Hamiltonians of 
the form 

H DG = A + kAx (1) 

where k is a coupling constant and Aq, A% are Onsager algebra fundamental generators. More generally, the 
Onsager (in)finite dimensional algebra in the basis of generators {A m ,G n }, m = 0, ±1, ±2, ... and n — 1,2,... 
reads 

[A n , A m ] = 4G m — n , [G m , A n ] = 2A n + m — lA n _ m , [G m , G n ] — . (2) 

In 5 , it was shown that Hjjg actually belongs to an (in)finitc family of mutually commuting conserved quantities 
provided the so-called Dolan-Grady relations 

[A , [A , [A , Ax]]] = 16[A , A x ] , [A x , [Ax, [A t , A }}] = 16[A U A ] (3) 
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are satisfied. As pointed out by Dolan and Grady, the power of the formulation (j2J relies in the operator statement. 
Indeed, Hamiltonians of apparently different quantum integrable models can be written as Q together with Q. 
For instance, it is the case for the Ising and XY models 5 , superintegrable chiral Potts models 6 and some 
generalizations [7j. 

Although the Onsager algebra (J2J) was the most important object in pp, it received less attention in the 
following years than the star triangle relations (which did not play any essential role in original Onsager's work) 
did. Actually, the most important progress in the approach of integrable systems was based on the star-triangle 
relations which originated in 018] and led to the Yang-Baxter equations, the theory of quantum groups, as well 
as the quantum inverse scattering method. This approach was further extended in [Sj: there, it was shown that 
the Yang-Baxter algebra has to be enlarged with the reflection equation (RE) which arises in various places: 
factorized scattering theory on the half- line |10j , integrable lattice models with non-periodic boundary conditions 
U3, non-commutative differential geometry on quantum groups. In this generalized framework, the integrability 
condition of a model follows from the existence of a solution to the reflection equation, the so-called reflection 
matrix. Indeed, it is the basic object which enters in the definition the quantum transfer matrix [5j denoted 
t(u) below. Using the fact that, by construction, [r(ui), r(u2)] = for all u\, ui and expanding in the spectral 
parameter u, the quantum transfer matrix provides a generating function for an (in)finite number of integrals of 
motion in the model. 

Whereas the inverse scattering method is well understood and provides many non-perturbative results in 
integrable models (spectrum, scattering amplitudes,...), a formulation based on an (in)finite dimensional algebra 
of the Onsager-type (J2J) would be obviously highly desirable: it would provide a new algebraic approach to 
integrable massive models 2 . Consequently, an interesting question is whether such algebraic structure might 
be actually hidden in the Yang-Baxter/reflection equations. One reason to believe so, is that similarly to the 
relations the explicit expressions of the R— matrix and if— matrix rely on the (usually non-local) quantum 
symmetry behind a model (lattice, quantum field theory,...) but not on its explicit microscopic (local) details. 
As all known models solved using the Onsager algebra's approach can also be analyzed in the framework of the 
inverse scattering method, it is expected that such a link can be exhibited. 

The purpose of this paper is to make a step towards the explicit construction of an (in)finite dimensional 
symmetry in massive integrable models, along the line initiated in |S]. Indeed, we will show that known 
integrable models with an underlying [7^1/2 (sfe) quantum group symmetry, for instance the sine-Gordon quantum 
field theory or XXZ spin chain, admit an alternative description 3 based on a "Hamiltonian" of the form 

H = A + A* (4) 

where q is a deformation parameter and the pair of operators A, A* (sometimes called a tridiagonal pair [26]) 
satisfies the "g— deformed" Dolan-Grady relations 

[A,[A,[A,A*y,-i] =p[A,A*] , [A*,[A*,[A*,A] 9 ],- 1 ]=p*[A*,A] . (5) 

Here we denote the q-commutator [A, B} q = q 1 / 2 AB — q~ 1 / 2 BA. Notice that these relations can be easily obtained 
from slightly more general ones called the tridiagonal relations (see also |14l IT5] ). Also, the special case q = 1 
in JSl leads to the Dolan-Grady relations J2J [I] (see also 00)- As we are going to see, the "Hamiltonian" 
differs from the usual local one, and is expressed in terms of nonlocal objects. Furthermore, it is the first of an 
infinite set of mutually commuting conserved quantities. On one hand, similarly to the undcformed case J2J) the 
relations © lead to an operator statement: they provide the integrability condition for models of the form 0} 
which does not depend on the dimension of the system or the nature of the space-time manifold. On the other 
hand, as we shall see the explicit realization/representation of the operators A, A* will depend on the model under 
consideration. 

2 The importance of such infinite dimensional symmetry might be compared to the Virasoro algebra with generators L n which play 
a crucial role in the analysis of conformal field theory 1111 . 

3 Except if explicitly specified, throughout the paper we assume generic values of q. Among the special cases, one finds the N = 2 
supersymmetric point in the sine-Gordon model which deserves special attention (for the case with a boundary, see 1121 1131 V 
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The paper is organized as follows. In Section 2, we derive a general solution (so-called K— matrix) of the 
RE in the spin— 1/2 representation of U q i/v{sl2). In a first part, assuming that the non-commuting entries 
Kij(u) of this A— matrix are Laurent polynomials of degree —2 < d < 2 in the spectral parameter u we obtain 
Kij{u) £ J-un(u; A, A*) where A, A* is a "Leonard" pair subject to q— deformed Dolan-Grady relations 10. More 
general if— matrices, denoted K^ n \u) below, are obtained using a standard "dressing" procedure. 

In Section 3, in order to derive in dependently the previous if— matrix we study the affine quantum group 
symmetry of the reflection equation. It is shown to be associated with a coaction tp acting on Leonard pairs which 
preserves the algebraic structure of the reflection equation and generates all K^- n \u) starting from K^°\u) = K(u). 
A set of relations for a (universal) if— matrix are also proposed. In some sense, they extend the well-known ones 
of Drinfeld for the R— matrix. For the spin— ^ representation and using this approach, we explicitly derive K(u). 
As expected, it agrees with the one obtained by solving directly the reflection equation. 

In Section 4, we consider the previous family of solutions K ( n ) (u) of the RE in the framework of the generalized 
inverse scattering method for the simplest case n = 0. We obtain naturally a generalization of the result of Dolan- 
Grady [S] for quantum integrable systems. Indeed, the corresponding Hamiltonian Q is a linear combination of 
A and A* satisfying (J5J, and enjoys a duality property exchanging A <-> A*. Its spectral problem is quasi-exactly 
solvable, except if q is a root of unity in which case it becomes exactly solvable. We give two simple explicit 
examples of representations for A, A* satisfying JSJ. For a single cyclic representation, the form Q coincides with 
the Azbel-Hofstadter Hamiltonian on a square lattice in the chiral gauge which spectrum is already known. For 
a single infinite dimensional representation, one finds Askey- Wilson polynomials as eigenfunctions of 

In Section 5, we give some explicit examples of two-dimensional quantum integrable field theory which contain 
non-local operators satisfying (JSJ). In a first part, we point out that some particular combinations of non-local 
conserved charges in the sine-Gordon quantum field theory actually satisfy the relations and in some special 
cases generate the Askey- Wilson algebra. It follows that an infinite set of non-local conserved quantities in 
involution can be obtained. The spectral problem for the simplest one J3J and higher conserved quantities is 
found to be quasi-exactly solvable and allows us to relate one-particle states to Askey- Wilson polynomials. 

In a second part, we apply more general cases n / of K^ n \u) to two-dimensional boundary quantum 
integrable models coupled with boundary degrees of freedom located at the boundary, which must satisfy JHJ) in 
order to ensure integrability. Related integrable models are, for instance, the boundary sine-Gordon model in 
the massive phase ^JEIE] or massless limit ^HIEO] coupled with various realizations of A, A*. It follows that 
the boundary space of state admits a basis in terms of Askey- Wilson polynomials. Finally, we briefly discuss the 
weak-strong coupling duality phenomena which relates the bulk algebraic structure (non-local conserved charges) 
to the boundary degrees of freedom. We find that both operators satisfy the same algebraic relations. Concluding 
remarks follow in the last Section. 

2 Quadratic algebras and deformed Dolan-Grady relations 

Among the known examples of quadratic algebraic structures, one finds the Yang-Baxter algebra. For further anal- 
ysis, let us first recall some known results. This algebra consists of a couple R(u), L(u) where the R— matrix solves 
the Yang-Baxter equation (YBE) and the L-operator satisfies the quadratic relations (with spectral parameters 
u,v) 

Rv V {u/v)L Vo (u) <g> L K (v) = L v;> (v) ® L Vo (u)R VoV;> [u/ v) (6) 

where Vb, Vq denote finite dimensional auxiliary space representations. Here, the entries of the L-operator act on 
a quantum space denoted V below. If one considers a two-dimensional (spin— i) representation for Vo and Vq, a 
solution of the YBE is provided by (in general, we will denote i?y v (u) = R(u)) 

R{ u ) = X] w y ( u ) G% ® a i ' ( 7 ) 

ije{o,3,±} 

where o~j are Pauli matrices, er± = (o~i ± ia2)/2 and are some combinations of elliptic functions. The 

corresponding L— operator in @ can be obtained in analogy of l|7)l. For inequivalent auxiliary and quantum 
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spaces, it reads 



L(u) = UJ i]i u ) a i ® Sj (8) 

J = {0,3,±} 

where the algebraic properties of the operators {Sj} are determined by the Yang-Baxter algebra JBJ. Here, 
the operators {Sj} act on the quantum space V. Comparing terms with different dependence on the spectral 
parameters u, v in © one arrives at a set of algebraic relations for {Si}, a quadratic algebra which is known as 
the Sklyanin algebra. However, a simpler solution of the YBE is obtained from the trigonometric limit of {u>ij(u)} 
in which case the R— matrix takes the form Q with 

u 00 (u) = I(gV2 + 1)(u _ g -i/2 u -i ) > ^{u) = \{q 1 ' 2 -l){u + q- 1 ' 2 u- 1 ) , 

"+-(«) = w-+(«) = <z 1/2 - <T 1/2 ■ (9) 

The corresponding limit of the Sklyanin algebra, also called trigonometric Sklyanin algebra (TSA), is also 
quadratic. However, it can be shown to coincide with the quantum enveloping algebra Ugi/z^sfa) with the identi- 
fication 

_ q S:i/2 + q- S3 / 2 q s 3 /2_ q -s 3 /2 

TTa I _1 / A 1 

where 



gl/4 + g-l/4 ' WJ g l/4_ q -l/4 



qs 3 _ ^-s 3 

[s3,S±] = ±5± , and 5 1 -] = /2 - /2 , (10) 

together with the Casimir operator 

w = + g-VSg-a + ( g l/2 _ q -Wf S -S + , (U) 

An other example of quadratic algebra is provided by the reflection equation (sometimes called boundary 
Yang-Baxter equation). This equation arises in various context (see for instance ^3 1^ UU G21 ) ■ Without spectral 
parameter, a systematic study of quadratic algebras associated with the RE can be found in |231 I24| . Here we are 
interested in the spectral parameter dependent form which appeared originally in |10| . For Vo = Vq it reads 

R(u/v) (K(u) ® J) R(uv) (J <8> K{v)) = ( J ® K{v)) R{uv) (K(u) ® I) R(u/v) . (12) 

Similarly to JHJ, in the spin— ^ representation Vo we introduce a if— matrix of the form 

K{u) = Oj- ®f2j(u) ■ (13) 

je{o,3,±} 

We assume that £lj(u) are Laurent polynomials of degree — 2 < d < 2 in the spectral parameter u with coefhcients 
in a (yet unknown) non-commutative associative algebra A. Replacing (|13(l in (|12(l . certain coefhcients are 
vanishing which leads to consider 

n (u) = (F + G)u/2- (G + F)u- l /2 , fi+(u) = Uu 2 + Vu- 2 + W , 

Q 3 (u) = (F-G)u/2- (G-F)u- 1 /2 , = Uu 2 + Vu' 2 + W (14) 

where {U, V, U, V, F, G, F, G, W, W} £ A. Replacing in (fT2"|) and after some manipulations we find for instance 

[17, V] = [u,w] = [v,w] = [u,v] = [U,W] = [V,W] = , 

[U,U] = [V,V] = , 

[U,F] = [U,G] = [V,F] = [V,G] = [U,F] = [U,G] = [V,F] = [V,G]=0, 

uv = vii, uv = vu , 
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and similarly substituting F <-> G, G <-> F. It immediately follows that {U, V, U, V} belong to the center of A. 
The fixed value of these elements over any quantum space representation V leads to introduce the linear relations 

U = c U , V = c V 

with co € (E. Then, from H12|l we also get 

(1 - q-^iUW - c UW) - [F, G] = , (q - l)(VW - c VW) - [F, G] = , 
FU - FV = FU - FV — GV — GU — GV - GU = , 

[F,G] -[F,G] = 0. (15) 

Assuming that W , are linearly independent, these relations imply 

V = q- 1 U 1 F = q~ 1 F G = qG. 

Let us now turn to bilinear expressions in F, G, W, W following from 112(1 . Using (|15p. they reduce to 

C/F(l - g ~ 2 ) + gVKG - GVK = , C/G(l - q 2 ) + - qFW = , 

c C/F(l - q~ 2 ) + gGM^ - WG = 0, c C/G(l - g 2 ) + FFF - ?fF = , 
[W,W] +(q-q- 1 )(q- 1 F 2 -qG 2 ) = 0. (16) 

The elements {W, W} e ^"un(F, G) may be any polynomials in {F, G}. However, coefficients of various terms are 
subjects to the remaining first equation of (|15() and eqs. (|16|l . Consistency implies that the allowed combinations 
are reduced to 

W = d[F,G] q -i +c 2 and W = -c a a[F, G] q + c c 2 , (17) 

where {ci, C2} G (F, together with the value of the element (in the center of .4) fixed to U = —\/c( j ciq~ 1 / 2 (q — q~ 1 ). 
It also follows that {F, G} are subject to the trilinear relations 

FG 2 + G 2 F -(q + q~ 1 )GFG - ^F = ^(q~ 1 / 2 - q l / 2 )G , 

C Cf Ci 

GF 2 + F 2 G -(q + q- v )FGF - -\g = -(q- 1/2 - q 1/2 )F . (18) 

COCJ Ci 

Notice that these relations are sufficient in order to satisfy the last eq. in i|ltj|) . Setting 

F = q 1/2 A and G = q- l ' 2 A* 

in (|18|l . it is interesting to observe that such kind of trilinear algebraic relations are called the Askey- Wilson 
relations first introduced by Zhedanov in |51| . They also appear in the context of representation theory associated 
with tridiagonal algebras © (for more details, see US]). According to an ordered pair A, A* of linear 
transformations on a finite dimensional vector space V 

A : V V and A* : V -> V (19) 

is a Leonard pair if there exists a basis for V with respect to which the matrix representing A is diagonal and 
the matrix representing A* is irreducible tridiagonal 5 , and similarly exchanging A «-» A*. A Leonard pair is 
essentially the same thing as a tridiagonal pair for which both operators are diagonalizable with all eigenspaces 

4 It is also possible to find infinite dimensional or cyclic representation V. Although examples can be found in 1261 and below, 
Leonard pairs in these cases remain to be define accordingly. 

5 A square matrix X is called tridiagonal whenever each nonzero entry lies on either the diagonal, subdiagonal or the superdiagonal. 
It is irreducible tridiagonal whenever each entry on the subdiagonal or superdiagonal is nonzero. 
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of dimension one [EH], Lemma 2.2]. The Leonard pairs are classified in [|25|. Theorem f.9] and Theorem 
5.16]. In it is shown that a Leonard pair A, A* satisfies 



A* A + A A* — (q + q~ )AA*A - j(AA* + A* A) - pA* = j*A + u>A + r\I , 
AA* 2 + A* 2 A-{q + q- 1 )A*AA*-j*{AA*+A*A)-p*A = 1 A* 2 +uoA*+r,*I (20) 

for a unique sequence of parameters 7, 7*, p, p* , oj, 77, if. To relate the quadratic algebra (|T2|) and ijT%|) to the 
algebraic relations lf2"U|) associated with a Leonard pair A, A* , we identify 6 7 = 7* = r) = rj* = and 

p = WS -^(q 1 /2_ g -l/2 ) and . . =Q (21) 

C Cf Ci 

where the parameters co,ci 7^ 0, C2 are arbitrary. Combining all previous expressions, we conclude that any 
solution of the reflection equation l|12|) of degree — 2 < d < 2 in the spectral parameter u - with non-commuting 
entries - can be written in the form (|13[) where 

fl {u) = (A + A*)(q l / 2 u - q~ l / 2 u~ 1 )/2 , 0, 3 (u) = {A - A*)(q l / 2 u + q~ l / 2 u~ l )/2 , 

a 1 / 2 ?/ 2 + q- x I 2 u- 2 
coCi(g-g x ) 

Q l / 2 U 2 + _1 / 2 M -2 

n_(u) = — — c ci[^, J 4*] 9 + coc 2 . (22) 

ci(?-« x ) 

More general solutions of the reflection equation l|12|) can now be obtained from i|13|) with l|22|l using the 
"dressing" procedure proposed by Sklyanin in [§] that we apply here: Knowing a solution of the reflection 
equation (|12J) . say K^(u), with quantum space V^ ) = V and given the Lax operator L j (u) given by (JHJ) with 
©, I|1U|I acting on a quantum space Vj we define a family of "dressed" reflection matrix 

K (n \u) = L n (uk) ■ ■ ■ Lxiuk^X^Lxiuk- 1 ) ■ ■ ■ L n (ufc _1 ) (23) 
acting on the quantum space 

n 

V (n) = (g) Vj <g> V . (24) 
j=l 

Provided [Sj,a] = for all a S A , it is easy to show that K^ n \u) satisfies the RE <|12|) for any value 7 of the 
parameter k. Choosing = K(u) defined by i|13[l with (|22() and (|20|l . i|21|) . we finally obtain a whole family of 
solutions of the RE JT2J. 

To conclude, let us mention that finite/infinite dimensional irreducible representations of K^> (u) in the spin— A 
auxiliary space representation Vo are classified according to the tensor product l|24|l . As the tensor product of 
representations Vj (associated with U q i/2 {si?)) is well understood, it is important to know all the Leonard pairs as 
well as their cyclic generalizations. In Section 4 and 5, we will provide explicit examples of such representations, 
related with integrable systems. 



3 Quantum affine reflection algebra 

Given any finite dimensional auxiliary space representation Vo = Vq (for instance the spin— \ representation) it 
is well-known that the (four dimensional) matrix solution of the Yang-Baxter equation 

R(u/v)R{u)R{v) = R(v)R{u)R(u/v) (25) 

6 Note that a Leonard pair, say A, A*, satisfying 121 H can be deduced from a different pair A, A* associated with 1201 for 7 = 7* = 
r] = rf = and using following substitutions in H20I : A — > A = A + al ,A* — > A* = A* + a*I and 7 — > — (q 1 / 2 — q^ 1 / 2 ) 2 a, 
p^ p + (q 1/2 - g- 1 / 2 ) 2 ^ 2 , uj ^ uj + 2{q 1 / 2 - q~ 1 / 2 ) 2 aa* , rj -* -pa* -uia- (q 1 / 2 - q- x l 2 ) 2 a 2 a* where {a, a*} G <E 

7 For k = 1 and up to an overall spectral parameter dependent factor one has L~^(u) ~ L(u _1 ). 
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can be derived using the underlying symmetry algebra. For the tensor product of two vector representations, the 
corresponding R— matrices were derived by Jimbo |29| for all non-exceptional quantum affine Lie algebras. In 
particular, using the U q i/2{sli) symmetry in the spin— | representation the R— matrix Q with JJjJ follows. 

An analogous method was initiated in (^Hj and further developed in 1321 ESj in order to derive some 
K— matrix solutions of the RE (|12|) using certain coideal subalgebras of quantum Kac-Moody algebras or Yangians 
32 . There arc, however, known examples of K— matrix derived using different symmetries that do not fit in this 
framework. For instance, in |18j we introduced an example of "dynamical" extension of coideal subalgebra of 
U q i/i{sl2). The purpose of this section is to reconsider from a more general point of view the meaning of this 
extension in the context of U q i/2 (SI2)— comodules. In particular, we propose a rather general approach to the 
construction of a universal K— matrix based on 1)40(1 . 



3.1 Preliminaries 

The universal R— matrix associated with the quantum Kac-Moody algebra [^1/2 (s^) was found by Drinfeld [3~4"| 
using the quantum double construction. This quantum Kac-Moody algebra, denoted for simplicity T below, is 
generated by the elements {hj, ej, fj}, j E {0, 1}. They satisfy the commutation relations 

qhi/2 _ q -hi/2 

[hi,hj] — 0, [hi,ej] = dijCj , [hi,fj] — —ciijfj , [ej, fj] — Sij ^ — ~~ Z1/2 
together with the q— Serre relations 

[e u [a, [e i} ej] g ] q -i] = , and [/*, [f h [/.,, fj] q ] q -i] = . (26) 

The sum k = ho + hi is the central element of the algebra. The Hopf algebra structure is ensured by the existence 
of a comultiplication A F ; T — > T % T and a counit £ j? : T — > (T with 

A^(e,) = e^g^ + ^/^e, 
A^(/i) - fi ® q h ^ 4 + q- h ^ 4 ® h , 

Ajr(hi) = hi®I + I®hi (27) 

and 

£r{ei) = £ F {fi) = £Ahi) = , E r {I) = 1 . (28) 
Furthermore, J- is a coalgebra [35] because 

(idjr X £jr) oAjr = idjr , (£jr X idjr) = idjr , (Ajr X idjr) O Ayr = (idjr X Ajf) O Ayr . 

Letting R E T ' ® T denotes the universal R— matrix, as shown in |34| it is the unique invertible element which 
satisfies the relations: 

A^{x) = RAjr{ x )R- x for all x E T , 

(id F x A F )(R) = R 13 R 12 , (A,; x id F )(R) = i?i 3 ^23 (29) 

where A'-p = a o A F and one introduces the permutation map a(x (g> y) = y ® x for all x,y E T . Here, using the 
definition R = ^2 i a i ®b l we denote R\2 = J2i a * ® & ® I, #13 = J2i I ®b l and i?23 — J2i A ® a % ® In 
order to obtain a solution of the Yang-Baxter equation (non-linear in the entries of i? for some finite dimensional 
representation), it is sufficient to find R E T ® T that satisfies (|29|l . Indeed, from the first equation of (|29|l one 
writes 

Ryi^T x idf){x ®y) = (A'jr x id?){x ® y)R\2 for all x,y E T . (30) 
Setting R = x <S) y and using the last relation of l|29|l one obtains the Yang-Baxter equation 
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R\2 R 13, R 23 = ^23^13^12 



(31) 



The universal version of © is derived similarly. One denotes the element L £ J- ® A, where A is a quadratic 
associative algebra. Actually, suppose that this quadratic algebra is such that there exists an element L = x ® a 
which satisfies (|30ll substituting idjr — > id*. Then L satisfies the Yang-Baxter algebra 

i?^2^1^2 = -^2^1^12 ■ (^2) 

Also, assuming a Hopf algebra structure for A with comultiplication : A — > .4 ® A, more general solutions 
lw of inm can be obtained. To construct them, one introduces the n— comultiplication 

A J n) : A — ► A &>■••<£> A , A l{ £ ] = (id A x • • • x id A x A^) o A J 1-1 * (33) 
for n > 3 with A^ = A^, and define 

L<») =(i^x A'|" ) )(i)=i n --- J L 1 (34) 

where indices here refer to the quantum space ,4. Acting with (i x A^ n) ) on (|32[1 and using the standard 

comultiplication rules, it is straightforward to show that L^> is indeed a solution of (|32|) . Obviously, we may 

i n \ 

have started from Aj^ instead. In that case, one shows that L = ■ ■ ■ L n also solves (13211 . In other words, 
the comultiplication A^ (or similarly A_^) preserves the algebraic structure associated with (13211 . 

Later on, we will focus on the spin— | representation of U q i/2 (sfo) and central element k — 0. Denoting 7ri 2 ' (x) 
the representation with spectral parameter u of the element x £ T, one has 

(tt^ 3 x Tr^ y )[R 12 } = R{u/v) , ] x id^L] = (35) 

and more generally 

(tt^ x id A x • • • x id$)[L (n) ] = L n (u) ■ ■ ■ L^u) . (36) 

In particular, the representations l|35|l are used in l|31|l and IMl'l) to obtain the Yang-Baxter equation (|25[l or the 
Yang-Baxter algebra ©, respectively. In the following subsection, we will proceed by analogy for the reflection 
equation. 

3.2 Comodule algebra and the reflection equation 

In this part, we study the algebraic structure related with the spectral parameter dependent reflection equation 
l|12l) . For F q (GL(2)) i.e. without spectral parameter, part of the following analysis was considered by Kulish and 
Sklyanin in [23] (see also [21], [121 and [37j). First, let us recall the concept of coaction and comodule necessary 
in the following analysis. Then, we propose a quite general framework in order to derive a "universal" reflection 
matrix. The special case of spin— i representation is treated in full details. 

Definition: |35j Given a coalgebra T with comultiplication A^ and counit Ef, A is called a left J- — comodule 
if there exists a left coaction map £ : A — > T <£> A such that 

(Aj, x id A ) o £ = (idjr x Q o ( , x id A ) o ( S id A . (37) 

Right T~ comodules are defined similarly. 

By analogy with the previous subsection, we are now looking for coaction maps if (and tp') which preserve the 
algebraic structure associated with the quadratic algebra {JT^J. Similarly to these coaction maps must be 

such that K^ n > is generated from K^ Q \ For a given (here spin— i) representation of J-, we previously obtained 
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the solutions K^ n \u) in terms of K^°\u). For instance, setting n = 1 in ((23|) and replacing k — ► v, it leads us to 
define the coaction maps by 

(7ri 2 ' x 7r„ x id^)[(«rfjr x <^)(-K")] = L(uv)K(u)L(uv~ 1 ) , 

(tt ( J } x tt v x id A )[{idjr x ip')(K)] = L(uv- l )K(u)L(uv) . (38) 

Here the notation tt v means that we do not specify the dimension of the representation for A. From now on, we 
identify A = T . Then, it is important to remember that the L— operators can be written in terms of the universal 
R— matrix as 

L(uv~ l ) = (ir$ ] x n v )[R 12 ] and L(uv) = {tt { J } x n v )[R 21 } (39) 
in order to introduce relations analogous 8 to 1(29|) : 

Conjecture: Let Ayr be the comultiplication associated with the coalgebra T . Let ip (and ip') be a two-sided 
coaction associated with an associative algebra A. There exists an element K £ J-®A which satisfies the relations: 

(i) <p'(a) = Kip^K^ 1 for all a e A, 

(ii) {Ayr x id A )(K) = (K ® I)R 21 (n ® K) , (A^ x id A ){K) = (M ® K)R 21 (K ® E) , 
{Hi) (idjr x <p)(K) = R 21 {K®I)R 12 , {idjr x <p')(K) = R\ 2 {K <g> M)R 2 \ . (40) 

Notice that the relation (ii) already appeared in [36 | . If such an element K exists, it is then straightforward to 
show that this element satisfies the reflection equation. There are two possible ways. For instance, from (i) we 
write 

(id? x (p')(x <g> a) (J ® K) = (1 ® K)(id?r x tp)(x ® a) for all x ® a e T ® A , 

or similarly from the first relation of (|29|l one has 

Rl 2 (A'jr x id A )(x ® a) = (Ayr x id A )(x ® a)Ri 2 for all x ® a G T ® A . 

Choosing K = x ® a and using (Hi) or (ii), respectively, in the two equations above we deduce the universal form 
of the reflection equation 

R 12 (K ® I)R 21 (I ® K) = (J <g> K)R 21 (K ® I)R 12 . (41) 

Similarly to the previous subsection, a generalized n— coaction ip( n > can be introduced for n > 2 with 
Together with l|38l) . it can be used to generate the family of solutions (|23|l starting from K(u) = K^°'(u). 
Consequently, the basic solution l|13l) with ll'l'li is a fundamental object. 

Here, we will not try to construct the universal reflection matrix. Instead, we will consider in details the 
derivation of K(u) in i|13|) in the spin— | representation of U q i/2(sl 2 ), using H4U|) . Note that a special choice of 
representation for the associative algebra A will not be necessary. The spectral parameter dependent form of R, 
K and L in the spin— ^ representation is given by (|35|l and 

(■k[? ] x ^ ] )[R 21 ] = R(uv) and (vri^ x id A )[K] = K(u) (42) 

from which (|f 2(1 follows. For higher dimensional representation one gets obviously the same form. 

Now, we would like to find the explicit expression of the action of ip on the elements of A. To do that, let 
us expand ((38(1 in terms of the spectral parameter v. There, we use the explicit form of K^(u) given by ((13(1 

8 However, a rigorous proof of the existence of a unique universal reflection matrix remains to be done. 
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together with JSJ. After some straightforward calculations, we find that the U q i/i{sl2) algebra arises explicitly 
through the evaluation homomorphism [/^^(s^) — ► U q i/i(sl2) 

7r„[ei] = vS+ , n v [e ] = vS- , 

Kv[fl] = V~ 1 S- , 7T„[/o] = iT 1 ^ , 

7r [ g fcl / 3 ]=d'» , n v [q h <>/ 2 } = q- s ° . 

The expansion of 1(38(1 finally gives 

(7ri 4) x n v x ^)[(^ x (^)(iv)]n = (tt^x/4 x ^)[gn 3 C(A)- g - 1 «" 3 C(^*) + C'Kw- 1 ;4,^*)] , 
(7rl |} x 7T, x «U)[(i<fr x ^')W] 2 2 = x i<U) [ gu 3 C( A* ) - g~ A) + C (u, w- 1 ; A* , A)] , 

(7r£ i) x tt, x x ip)(K)] 12 = (ir vql/i x id A )[- ± g "'"^ (00 + O V~ 2 ; A A*)] , 

coc x (g-g x ) 

(tt!' 3 x tt„ x t7f,i)[(i(fr x (^)(X)] 2 i = K,!/4 x id A )[- qU ( +q _ ClO + C'V^'V*,^)] (43) 

where 

C(J) = J ® / , 

C(A) = (c+ei + C-h)q h ^ A ® J + q' ll/2 ® A , 

C(A*) = (c_e + c+Zo)?' 1 "/ 4 ® 7+ q ho/2 ® A* . (44) 

The terms C"i C do not contribute in the asymptotic u — * oo or u — > so we don't write them explicitly 
here. The parameters c+, c_ read 

c+ = -l/coc^g 1 / 2 + g- 1 / 2 ) and c_ = -l/c^g 1 / 2 + q' 1 ' 2 ) . (45) 

It should be stressed that each of these maps £, £" define a left coaction and confirm the fact that A is a left 
U q i/2(sl 2 )~ comodule. For instance, it is easy to check using (J27J), (|25)l that ( defined by is consistent with the 
comultiplication A^- = a o Ajf, i.e. satisfies (13711 . Starting from (i) in (J3DJ, we can now write (the representation 
in v— space is not specified yet) 

<-) (-) 
(iTu x n v x id A )[{E ® K)(idyr x ip)(x ® a)] = (7r„ 2; x 7r„ x id A )[(idjr x y> )(a: (8 a) (J ® if)] (46) 

for all x G JF, a G A. Setting K = x ® a and using (m) of (|4t)(l together with 1)39(1 . 1(38(1 . both sides can be 
expanded in the spectral parameter u. In the spin— ~ representation, setting u — ► oo or u — > it follows that the 
leading terms are given in terms of ((44(1 . We find that the if— matrix 1(13(1 must satisfy 

K{v)n%M*)]=*% q y»)]K{v) for all a G {I, A, A*} . (47) 

This intertwining equation is actually sufficient to determine K(u). Indeed, in this representation one has 
7r^)(5 ± ) = a± , 7P')(s3) = 03/2. From the diagonal part of 1(47(1. it is easy to get the expressions flo(v) 
and ^(u) in ((22(1 . Replacing them in the off-diagonal part of 1(4711 . it follows that Kij(v) G Funiv; A, A*) for 
i 7^ j are polynomials of degree —2,2,0 in the spectral parameter v. Also, one finds easily that the spectral 
parameter dependent terms must commute with A, A*, whereas the remaining terms (denoted W,W in 1(14(1 ) 
must satisfy for instance: 

q l ' 2 A*W - q-^WA* = -C-Aiq 1 ' 2 +q- 1 ' 2 ) , 
g-^AW" - g 1/2 fi = c_A*(g x / 2 + g" 1/2 ) (48) 

and similarly for W using c+. However, from the diagonal part of 147(1 . one also finds c+(W — CqW) = [A, A*]. 
Replacing in 1(48(1 one immediately gets the trilinear algebraic relations 1(20(1 together with ((21(1 . Consequently, the 
solution in the spin— i representation 1(13(1 with 1(22(1 . 1(20(1 . ((21(1 can be derived directly from (i) in 1(40(1 as soon 
as one knows the explicit form of the coaction acting on A, or at least its asymptotic behavior in the "auxiliary" 
spectral parameter. To resume, the coaction plays here a role analogous to the coproduct in case of a Hopf algebra 
structure. 
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4 General features of the integrable structure 



Either considering the constraints which follow from the reflection equation, or studying the quantum affine 
reflection symmetry, we have seen that the trilinear algebraic relations (|20|l with (|21|l arise in both approaches. 
However, it is clear that the parameters Co, ci, ci entering in the reflection matrix 1)13(1 are not restricted by any 
symmetry argument during this process. Instead, their explicit values will depend on the explicit realization of 
the Leonard pair A, A* i.e. on model- dependent and representation characteristics of the corresponding operators. 
These realizations, as we will see below, admit finite/infinite dimensional or cyclic irreducible representations. 
Consequently, to speak in full generality i.e without specifying a realization/representation one should instead 
of 120() consider the q— deformed Dolan-Grady relations © as the integrability condition for the class of models 
associated with (fl2*|) . 

Having identified this integrability condition, the next step is the construction of related integrable models. 
In the quantum inverse scattering framework, there are many examples of quantum integrable systems which 
are constructed starting solely from a L— operator solution of the Yang-Baxter algebra ©. Using the coproduct 
structure, it follows that the entries of (|34() (see for instance in the Appendix) are rather nonlocal objects, usually 
expressed in terms of the local degrees of freedom of the system acting on the quantum space. Taking the trace 
over the auxiliary space of the monodromy matrix, one obtains among conserved quantities the Hamiltonian of 
the system. By comparison between L and K i.e. JBJ and 1)121) . we can then consider a quantum integrable model 
constructed solely from a K— operator of the form (|13fl with 1)221) i.e. choosing n = in l)23|) . Its generating 
function is given by 

T(u)=Tr Vo [K(u)] . (49) 

Explicitly one finds t(u) — [uq 1 / 2 — u~ 1 q~ 1 / 2 )Y\ where the Hamiltonian 9 takes the form (@J with A — > A : 
A* — > A* which satisfy the q— deformed Dolan-Grady relations JHJ. Beyond the Hamitonian, note that there 
exists a Casimir operator 

Q = qp* A 2 + q- 1 P A* 2 - 3- [A,A*] q [A,A*] r i + ^ g ~ g — u>(AA* + A* A) . (50) 

An interesting problem is now to construct a family of quantum integrable models associated with (0J in 
the spirit of |S] based on more general objects, say A, A*, satisfying JSJ but not l|2U|l . First, by analogy with 
the undeformed case @] it is indeed well expected that one can find operators A, A* which admit two types of 
representations: 

(I) - "Single" representations i.e. A, A* are represented by a pair A, A* which obey l)2UJI: 

(II) - "Multiple" representations (non-trivial tensor products of irreducible ones) i.e. A, A* satisfy JSJl, not (|20|l . 
Note that representations of type (II) can be obtained from more general solutions 1)23)1 of the reflection equation, 
as shown in 001- Secondly, the power of the algebraic constraint J5J) relies on the fact that it is an operator 
statement which does not refer to any local structure of the model under consideration. It is thus natural to 
expect that A, A* are rather nonlocal objects. By analogy with the undeformed case jo], it is possible to show 
that there exists an (in) finite family of mutually commuting conserved quantities of the form (r > 0) 

G (2r) (A,A*) = A (2r) +A* (2r) with G (0) (A,A*) = H (51) 

provided © is satisfied. For example, assuming 10 it is not difficult to check that 

G ( 2 )(A,A*)= [A > [A l A*] g _ l ] g +[A*,[A*,A] fl _ 1 ], (52) 

is conserved. We report the detailed construction of higher conserved quantities to a separate work [40) . Then, an 
important problem concerns the diagonalization of these conserved quantities G^ 2r \A, A*). It is actually closely 
related with the problem of partial algebraization of the spectrum associated with difference equations which 
arises in the context of quasi-exactly solvable systems. Two situations arise: 

9 Obviously, using a scale transformation of Leonard pairs one can introduce extra parameters related with p, p* in CT 
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(1) q is root of unity: the spectrum associated with G^ 2r '(A, A*) is finite. All periodic eigenfuntions can be 
obtained algebraically; 

(2) q is not a root of unity: the spectral problem associated with G( 2r )(A,A*) is quasi-exactly solvable. It 
is related with the representation theory of (JSJ). Only a part of the spectrum can be obtained, with polynomial 
eigenfunctions (for instance Askey- Wilson, big g-Jacobi polynomials,...) or their restrictions |26| . 

Using the explicit realization of the operators A, A* in terms of known algebras (see below), for single rep- 
resentations (type (I) above) the spectral problem associated with (|5T)l for r > can be reduced to a single 
^—difference equation of the form 

a(z)$(qz) + d{z)^(q- 1 z) - v(z)$(z) = Atf(z) (53) 

where a{z), d(z), v(z) are polynomials in z defined according to the model under consideration, the realiza- 
tion/representation of the Leonard pair and the explicit form of G' '. The eigenfunctions 9(z) of the algebraized 
part of the spectrum are also polynomials of finite order, Askey- Wilson polynomials for instance. Then, A denotes 
the corresponding subset of eigenvalues of the operator G' ' (A, A*). 

To be more explicit, let us give two examples corresponding to the cases (1) and (2) above and single rep- 
resentations (type (I)). For instance, following our results in ^S] it is not difficult to find a realization for the 
operators A, A* in terms of the Weyl algebra. Choosing 

1/2 *l/2 

A = -*-£ riQ + Q' 1 ) and A* = -t-£ -(P + P' 1 ) where PQ = q~ 1 QP, (54) 

q-q 1 q-q 1 

one can check that A, A* satisfy J^J and also 1(20(1 for 7 = 7* = lu — 77 = 77* = 0. For q = exp(2i%M/N) 
root of unity where M, N are mutually prime integers, it is well-known that the Weyl algebra can be realized 
using N x N matrices with the matrix elements {m\Q\n} = q~ m 8 m ^ n and (m\P\n) — <5 m +i. n [mod N) f° r which 
one has P N = Q N = 1 . In this case, the pair A, A* admits a finite dimensional (cyclic) representation but the 
matrix representing A* is not tridiagonal. For these values of q, the Hamiltonian (@J together with 154(1 actually 
coincides with the Azbel-Hofstadter (AH) Hamiltonian 02j which describes the problem of Bloch electrons in a 
magnetic field on a two-dimensional lattice. This model has been studied in details in several works 10 (see for 
instance (431 1441 1431 146| ) . Using the relation between the group of magnetic translations and the quantum group 
Uq{sl2), it was shown that the AH Hamiltonian can be written as a linear combination of the quantum group 
generators |43| . Also, its spectrum can be represented in terms of solutions of Bethc ansatz equations on high 
genus algebraic curves (431 I44| . In particular, for q a root of unity the polynomial eigenfunctions and the Bethc 
ansatz cover all the spectrum. Previous analysis shows that the AH Hamiltonian on a square lattice in the chiral 
gauge can be written as a linear combination of operators of the form 1)54(1 ■ the operators P, Q being identified 
with the generators of magnetic translations in each direction of the two-dimensional lattice. Using the bilinear 
formulation of ((20(1 in terms of the Askey- Wilson algebra [38] , it also explains the relation between the zero mode 
wave function and AW polynomials pointed out in |43| . Consequently, the AH Hamiltonian provides an example 
of quantum integrable system associated with the q— deformed Dolan-Grady relations (JSJ. 

Our second example is quasi-exactly solvable i.e. corresponds to generic values of q. For simplicity, let us 
consider the values p = —(q — q^ 1 ) 2 and p* = —abcdq^ 1 (q — q~ r ) 2 in J^J. Leonard pairs in this case admit an 
infinite dimensional representation in the basis of the Askey- Wilson polynomials with continuous weight 11 |2fij . 
These q— orthogonal polynomials, here denoted P n (z) with integer n > 0, are symmetric Laurent polynomials of 
the variable z i.e. they are invariant under the change z — ► z _1 . They can be written explicitly in terms of the 
basic q— hypergeometric functions as 

„ , , _ / q~ n , abcdq™" 1 , az, az' 1 \ . . 

P ™^ = 4 M at , ac ad " ( 55 ) 



1(> For an alternative derivation of the AH Hamitonian starting from a L— operator instead of K, we report the reader to the 
Appendix. 

11 For finite dimensional representation, it leads to Askey- Wilson polynomials with discrete weights. 
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In particular, these polynomials are known to satisfy the three terms recurrence relation 



p n P n+ i(z) + a n P n (z) + 7 „P„_i(z) = [z + z r )P n {z) 



(56) 



where P_! = 0, Pi = 1 and 



In 



(1 - abq n )(l - acq n )(l - adq n ){\ - abcdq 11 ' 1 ) 
a(l - abcdq 2n - r )(l - abcdq 2n ) ' 

a(l - q n )(l - bcq™- 1 )^ - bdq n - r )(l - cdq n ~ x ) 
(1 - abcdq 2n - 2 )(l - abcdq 2 ^ 1 ) 



a + a - (3 n - In ■ 
In this basis and choosing the symmetric case abed = q, one has 

\ 



tt(A) 



/ a 7i 

Po "1 72 

Pi a 2 

V 



7T(A*) =diag (2,q + q- i ,q 2 + q- 2 ,-,-, 



(57) 



On the other hand, these polynomials are eigenvectors [HHI of a second order q— difference operator ID with 
parameters a, b, c, d: E)P n (z) = (abcdq n ~ 1 + q~ n )P n (z). The Leonard pair in this representation becomes 7r(A) = 
z + A* = ID. Setting abed = q, it follows that the Askey- Wilson polynomials (|55J) satisfy the q— difference 



equation 



where 



ttz)P n {qz) + az-^PniQ- 1 *) = (£(*) + - (1 - ?"")(! - g"))P n («) 



(58) 



(1 - az)(l - bz)(l - cz)(l - dz) 
(1 - z 2 )(l - qz 2 ) 



Using (|5(j|l . I|57|l and 1)58(1 . the spectral problem associated with the Hamiltonian Q for p = p* 
takes the form 1(53(1 together with 

a(z)=(i(z), diz)^^' 1 ) and u(z) = £(z) + C(z _1 ) - z - z" 1 - 2 . 



(59) 



Using the standard Bethe ansatz technique, it follows that the roots zi"' of the Askey- Wilson polynomial l(55|) - 
eigenfunctions of 10} - determine the spectrum and obey the Bethe ansatz equations. This problem was treated 
in details in so we don't need to repeat the analysis here and report the reader to this work for details. 

Finally, it should be stressed that according to the model and realization of the Leonard pair considered, the 
values of the structure constants p, p* , 77, r\* , u> in 1(20(1 determine the form of polynomial eigenfunctions. Obviously, 
the two examples above do not exhaust all possibilities. For instance, if u> ^ and p = p* = 0, the eigenfunctions 
reduce to q— Hahn polynomials associated with 3<E>2- For uj = p = p* = they correspond to q— Krawtchouk 
polynomials associated with 2 < E > i ■ 



5 Related two-dimensional quantum integrable models 
5.1 The sine-Gordon model revisited 

For the coupling constant p 2 < 2, the bulk sine-Gordon model in 1+1 Minkowski space-time is massive and 
integrable. The particle spectrum consists of a soliton/antisoliton pair (tjj + (9),ip-(8)) with mass M and neutral 
particles, called "breathers", B n {ff) n = 1,2,...,< A. As usual, E = M cosh6* and P — Msinhfl the on-shell 
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energy and momentum, respectively, of the soliton/antisoliton. In this model, TL is the Fock space of multi-particle 
states. A general TV-particles state is generated by the "particle creation operators" A ai (9i) 

\A ai (6 1 )...A aN (6 N )) = A ai (6 1 )...A aN (6 N )\0) (60) 

where characterizes the type of particle. The commutation relations between these operators are determined 
by the S-matrix elements, which describe the factorized scattering theory |47|. Integrability imposes strong 
constraints on the system which imply that the general S'-matrix factorizes in two-particle/two-particle amplitudes 
which satisfy the quantum Yang-Baxter equations. To determine this basic S'-matrix, an alternative approach 
was proposed in 0S| based on the existence of non-local conserved charges in quantum integrable field theory. For 
the sine-Gordon model, the symmetry algebra is known to be identified with U qo (sl2) where the relation between 
the deformation parameter and the coupling constant is go — exp(— 2iir/(3 2 ). To exhibit this underlying quantum 
symmetry, the authors of |48) used the Lagrangian representation of the sine-Gordon model: 



Asa = I d 2 x 



(^(5^) 2 -2 Ai cos(/30)) . (61) 



Here ijl is a mass scale and <fi is the sine-Gordon field. In the perturbed conformal field theory (CFT) approach, the 
corresponding free bosonic fundamental field is decomposed in its holomorphic/antiholomorphic components 12 

4>(x,y) = ip(z)+Tp(z) . 

The non-local conserved currents J±{x,t), J±(x,t) of the form 

2i - 2i 

J±{x,t) =: exp(±— <p(z)) : and J±(x,t) —: exp(^— ^(z)) : 
(3 (3 

generate the non-local conserved charges denoted Q±, Q± below. Due to the non-locality of the currents, at equal 
time t they possess non-trivial braiding betweem themselves. Then, together with the "bulk" topological charge 

T = L [ dxd x cb (62) 

they generate the quantum enveloping algebra U qo (sl2) |48j. Contrary to the local integrals of motion, these 
conserved charges are not in involution. Their action over asymptotic soliton states can be deduced from a field 
theoretic approach for a single asymptotic state 13 with topological charge 2m 

\v;m) where v = exp(0) , — j < m < j (63) 

in the spin— j representation and fixed rapidity 9 it reads 

tt„[Q±] =ce xe S±q± S3 , tt„(Q ± ] =ce- xe S±qT S3 and n v [T] = 2s 3 (64) 



where 

A = 4- - 1 and c A = -i-^(qg Z - 1) 



— 1 and c 2 = — i-rr;(" 

/3 2 A 2 



More generally, the action of the non-local conserved charges Q±, Q± and T over a multi-particle state is obtained 
using the coproduct rules of U qo {sl2) (see @H1 for details). In the following, we will associate one-particle states 
to single representations (type (I)) and multi-particle states to representations of type (II). 

12 If wc denote the expectation value over the Fock space vacuum of the CFT by {...}, then these components are normalized such 
that (ip(z)ip(w)) = — ln(z — ui), ((p(z)(p(w)) = — ln(S — w), (ip(z)ip(w)) = 
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For the fundamental soliton/antisoliton with topological charge ±1, j = 1/2. Higher values of j occur for boundstates. 
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5.1.1 Askey- Wilson dynamical symmetry and duality 

A natural question is wether certain combinations of non-local conserved charges might satisfy the relations JSJ 
and generate conserved quantities in involution. In order to exhibit such non-trivial combinations in the sine- 
Gordon model, we use the following trick. According to the analysis of Section 2 and (|23[) . the simplest (up to 
an overall factor) non-trivial "dynamical" solution of the RE (|12l) is of the form (|13fl with l|22l) . An easy way to 
find single representations (type (I)) of operators satisfying (|SJl is then to start from a non-dynamical solution 
K(°\u), even very simple, and to apply (|23|l for n = 1 using the Lax operator For instance, let us choose the 
simple off-diagonal constant solution of (|12|) given by K = (c + cr + + c_a_)/(g 1 / 2 — q^ 1 ^ 2 ) and set v — > wo9 -1 ^ 4 hi 
<|38|) . It is straightforward to obtain the solution with (|22|l through the identification 

A = c + v S+q S3 / 2 + c_v- 1 S^q S3 / 2 and A* = c_v S-q~ S3 / 2 + c+v^ 1 S + q~ S3 / 2 . (65) 

Here we use l|45[) . <|20|) with <|21[1 and we find the relation 

where iUj is the Casimir operator (jll[) eigenvalue 

^ = q J+1 ' 2 + g^'- 1 / 2 . (67) 

An other useful realization of operators satisfying (J5J can be obtained from the non-diagonal non-dynamical 
solution of the reflection equation 14 [23| , [113 • It yields to 

A ee c+^oS+g" 372 + c-V^S-if afi + €+q 33 and A* = c_u S'-<r S3/2 + c+v^S+q-" 3 ' 2 + e-q-* 3 (68) 

where e± are free parameters. Actually, it is not difficult to check that the above simple realizations of A, A* 
satisfy (|20|l with 7 = 7* = 0, p = p* = 1/cqc 2 and the structure constants 



c 2 (^o; j) + e+e-CiQ? 1 / 2 - q 1 / 2 ) 1/2 _ -i/2\ 

Cl 



e+ KV 1/2 + ^o 2 9 1/2 ) 



V*( v 0;j) = V( v 0'J)\e+^e- ■ (69) 

Here, en, ci are non- vanishing parameters, C2 is the same as above ()66|) and wj is the Casimir eigenvalue (jllfl . It 
is now straightforward to identify the operators satisfying (0 in the sine-Gordon model. First, setting 

v Q = exp(A6») , q 2 = q , c + c = c_ ; 



d(gV2 + g -i/2) 

and using H64|l one can interpret H68|) as combinations of non-local charges acting on a single (one-particle state) 
representation H63|) i.e. of type (I). In particular, the fact that l|68[) satisfy the trilinear algebraic relations (|2UH 
with l|69[) implies ©. More generally, in order to find A, A* which do not reduce to a Leonard pair, one has 
to consider the action of the non-local conserved charges over multi-particle states (type (II)). In this case, we 
propose to consider the pair 

Q+ = —Q+ + Q- + e+qZ and Q_ = Q_ + — Q+ + e-q T . (70) 
c c 

To show that q— deformed Dolan-Grady relations are satisfied for any N— particles state representation of l|7U|l it 
is necessary and sufficient to show that the coproduct of U qo (sl2) leaves the relations (JSJ unchanged. Thus, we 

For a similar calculation, see \M . Note that this non-dynamical solution can also be obtained from 1221 for A, A* constants. 
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have checked explicitly that the non-local conserved charges satisfy the q— deformed Dolan-Grady relations JSJ 
with the substitution 

A -> Q+ and A* -> Q_ (71) 

for any Co 7^ and free parameters e± (for details, see @U|). I n this case, the value of p, p* in JSJ is 

. (l + g )(l + g-l) 

p = p = . (72) 

co 

We would like to stress that the trilinear algebraic relations H2()(l only arise in case of a single representation 
(one-particle state). This algebra actually admits a bilinear formulation in terms of the Askey- Wilson algebra 
AW(3) introduced and studied by Zhedanov in [HI]. It is generated by the elements {Kq, K%, K2} which satisfy 

[K , K x ] q = K 2 , [K 2 , K ] q = BK Q + C X K X + D 1 I , [K X ,K 2 ] q = BK X + C K + D I 

where B, Co, Ci, Do, D x are the structure constants. It is an exercise to relate these structure constants to the 
ones in l|t)9[) . According to the analysis of |51| . one can obtain irreducible finite dimensional representations 
depending on the value of the rapidity and the dimension of the asymptotic state representation 2j + 1. Following 
|51| . it is possible to find a basis \i/j m ), — j < m < 3, f° r the asymptotic particle state of rapidity 9 which 
diagonalizes tt v [Q+] ■ On the other hand, one can find a basis \<p s ), —j < s < j for the asymptotic particle 
state which instead diagonalizes tt v [Q-\ ■ Consequently, the conserved quantities (|70|l enjoy a remarkable duality 
property which reminds the one pointed out by Dolan-Grady in their paper [S]. Since both representations have 
the same dimension, both basis are related by a linear transformation \ip s ) — ^2{s\p) \ip p ). As shown in |51j . the 
corresponding overlap functions (s\p) can be explicitly expressed in terms of the Askey- Wilson polynomials on a 
finite interval of the real axis. 



5.1.2 Related spectral problem and asymptotic one-particle states 

It is now possible to construct an infinite set of mutually commuting conserved quantities of the form (|51(l in the 
sine-Gordon model (for a detailed analysis of a slightly more general integrable structure, see 001); the first ones 
being 

H = Q+ + Q-, 

G (2) - [Q + ,[Q + ,Q-] q _ 1 ] q +[Q-,[Q-,Q + ] q _ x ] q . (73) 

In general, the corresponding spectral problem is rather complicated. However, acting on a one-particle state 
all expressions drastically simplify. The corresponding spectral problem for a state of rapidity 9 and topological 
charge 2m in a representation of dimension 2j + 1 reads for r = 0, 1,... 

G^ 2r ^ \v; m) = A2r(v; j, m)\v; m) for all — j < m < j . (74) 

Up to an overall (v, j)— dependent coefficient which changes according to the value of r, this spectral problem 
always reduces to the one for H. Indeed, using the trilinear algebraic relations l|20(l it is not difficult to check 
that if \v; m) is an eigenstate of H, so it will be for all higher conserved quantities. For general values of qo, it is 
a quasi-exactly solvable problem. To study it, we follow the analysis of 02 . Expressed in the weight basis, the 
generators of U qo (s^) leave invariant the linear space (of 2j + 1 dimension) of polynomials F(z) of degree 2j. The 
lowest/highest weights are identified with Fq = 1 and F 2 j = z 2j , respectively. One has 

qt'F(z) = q^F(q^z), 

S+F(z) = - Z (qo 2j F(q z) - ql j F{q^z)) , 
(?o -Qo ) 

S-F(z) = 1 (F(q a z)-F(q^z)) . 

z(qo - q ) 
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In this basis of polynomials, for general values 15 of the parameters we may represent a one-particle state with 
rapidity 9 (recall that v — exp(0)) by 

2j 

(z\v;m)=V m (z) with * m («) = JJ (z - z™) (75) 

n=i 

where zi™' denote the roots of the m— th polynomial. From l|74|) and using the representation above, it is 
straightforward to obtain a q— difference equation of the form l|53|l with the substitution q — * q2o where the 
coefficients are explicitly given by 

a(z) = (-c+vzq^ + C-V^ 1 z^q^ 3 + e + (q - q Q 1 )q^ 23 ) /(q - q^ 1 ) , 
d(z) = ( - c_vz~ 1 q 3 ) + c+v^zql 3 + e_(q - q^q^) / (q - q^ 1 ) , 

v(z) = (- c+vzq> Q + C-V^z^q^ 3 + c + v~ 1 zq G 3 - C-VZ~ 1 q 3 Q ) / (q - q^ 1 ) . (76) 

The corresponding eigenvalue Ao(v;j, m) can be deduced from the results of |41|. It is expressed in terms of the 
roots z^ in (|75|1 which have to satisfy a system of 2j algebraic (Bethe) equations. This system has exactly 2j + 1 
solutions corresponding to 2j + 1 eigenfunctions ^> m (z), —j<m<j. Note that these roots depend on value of 
the rapidity 9 and e± . 

To resume, the sine-Gordon model possesses an underlying dynamical symmetry associated with (JjJJ with 
H71[). For single representations, it is isomorphic to the quadratic Askey- Wilson algebra, which is generated by 
the charges. This symmetry is preserved at all orders in perturbed CFT framework because the fundamental 
non-local charges Q±,Q± and T are conserved. The point is that these charges generate an infinite number 
of mutually commuting conserved quantities G^ 2 ^ |4U| . Actually, their actions on asymptotic states mix non- 
trivially the holomorphic and antiholomorphic sectors (in the perturbed CFT approach). Then, conservation of 
G^ for any r imposes a severe restriction on the asymptotic states structure. For instance, a one-particle state 
l|75|l admits a representation in terms of polynomials whose roots satisfy Bethe equations. The corresponding 
eigenvalues can be computed using the usual Bethe ansatz techniques (see for instance |41|). 

5.2 The sine-Gordon model with/without a dynamical boundary 

Following 9 , boundary integrable models on the lattice - for instance the XXZ spin chain - can be obtained using 
l|23|) . In ^7] we have applied Sklyanin's formalism [3] to construct a classical integrable field theory restricted 
to the half-line, namely the sine-Gordon model, coupled with a mechanical system at the boundary. In general, 
such two-dimensional classical field theory will remain integrable provided there exists a classical solution K c \ (u) 
of the classical reflection equation. In this case, the generating function reads 

t c1 (u) = Tr Vo [T^K^f^u' 1 )} (77) 

where T c i{u) (and similarly T c i(u)) is the sine-Gordon classical monodromy matrix. Expanding in the spectral 
parameter u one obtains an infinite number of mutually "Poisson" commuting quantities, which give the integrals 
of motion and ensure the integrability of the classical system with a boundary. However, it is clear that the analysis 
of J7| for the classical boundary sine-Gordon model can be extended to more general boundary operators provided 
they satisfy some classical quadratic algebra. This can be done starting from the classical limit of the reflection 
matrix Expanding t c i(u) in the spectral parameter u, among the integrals of motion one identifies the 

Hamiltonian. Its boundary contribution is found to be linear in the classical boundary operators coupled with 
the sine-Gordon classical field 4>. Setting the coupling constant of the sine-Gordon model $ — 0, this boundary 
term reduces to the classical analogue of 

At quantum level, an integrable (massive) boundary sine-Gordon model coupled with a quantum mechanical 
system at the boundary was proposed in |18j . There, the asymptotic behavior of the boundary operators was 

15 For some special values of the parameters, the degree of the polynomial eigenfuntions can be different. Here we do not consider 
this possibility. 
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expressed in terms of the operators (|54|) . However, as we have shown in Sections 2 and 3, the existence of a 
K— matrix is related with i|2U|) (see also 18 for details). It is thus more general to consider in full generality 
boundary operators which asymptotically satisfy l|2l)|) rather than choosing a (somehow simpler) realization like 
(J3J. Then, we propose 

A bsg =J dy J dxQ-(d v 4>) 2 -2^cos(/30)) + V 1/2 J dy $* ert (y) + A boundary (78) 

where the interaction between the sine-Gordon field and the boundary quantum operators reads 

^ ert (y) = £_(yy^y^ + £ + (y)e~^ ^ 

and Aboundary is the kinetic part 16 associated with the boundary operators £±(y). In the deep UV, this model 
can be considered as a relevant perturbation of a conformal field theory (for instance the free Gaussian field) on 
the semi-infinite plane with dynamical boundary conditions at x = 0. In this case, similarly to the bulk (|62(l the 
free Gaussian field restricted to the half- line can be also written in terms of its holomorphic/anti-holomorphic 
components. However, these components are instead normalized such that 

(cp(z)ip(w)) = -21n(z - w), (ip(z)ip(w))o = -21n(z - w), (ip(z)(p(w)) = -21n(z - w) 

where (...)o now denotes the expectation value in the boundary conformal field theory (BCFT) with Neumann 
boundary conditions. The boundary conditions in Ij78l) can be derived from the expectation value of the local field 
d x (t>(0, y) with any other local field in first order of boundary conformal perturbation theory fi — > 0. Actually, 
these quantum boundary conditions take a form similar to the classical ones |18| . 

We can choose the y-direction to be the "time" in which case the Hamiltonian contains the boundary con- 
tribution, and the Hilbert space Hb is identified with the half-line y = Const.. Then, correlation functions are 
calculated over the boundary ground state denoted |0)b below. This state can be expanded as 

|0) B = \vac) B <8 \0)bcft + 0{fi) , (79) 

where \0)bcft <= "Hbcft belongs to the Hilbert space of the BCFT and \vac)s is an eigenstate of the boundary 
Hamiltonian for the coupling /3 — * 0. 

There are good reason to believe ^3 E| t na t, a quantum analogue of the classical integrals of motion can be 
constructed explicitly. Similarly to the bulk case, instead one can show that the model l|78() remains integrable at 
quantum level too using the existence of boundary non-local conserved charges. Along the line of |3()l 1321 IT%] it is 
possible to construct non-local conserved charges in (|78(l in terms of the bulk ones (which are no longer conserved 
individually). Assuming the same asymptotic behavior at the end point of the time axis 

S + (y = ±oo) ~ £+ and = ±oo) ~ £- , (80) 

it follows that the quantum affine reflection symmetry of the model l|78|l is generated by the operators 

0i 0) = Q± + Q T + £±1o Tb with r b = I dx d x <f> (81) 

and £±{y) = fj}^ 2 £±(y)$ 2 /(l — p 2 ) ■ For the special case (I54f) . this was shown in details in Here, % denotes 
the boundary topological charge which is no longer independently conserved. In the following, we are going to 
consider different restrictions of the boundary operators. This will reveal interesting features of the model l|78l) . 
In particular, most of the known and well-studied massive or massless boundary integrable models can be seen as 
limiting cases of (JTSJ- 

• The boundary sine-Gordon model revisited. The quantum boundary sine-Gordon model introduced 
by Ghoshal-Zamolodchikov is obtained for the one-dimensional (trivial) representation of the boundary operators 

16 A solution of the "dual" reflection (associated with a boundary on the left hand side) can be used to introduce such term, without 
breaking integrability. 
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(which become free parameters) in (|78f) . In other words, the boundary degrees of freedom are "frozen" to 
constants. The corresponding non-local conserved charges calculated in (31)1 132] take the form (|81|l together with 
the substitution £± — > e±, where e± are free parameters 17 . It is then interesting to notice that, from an algebraic 
point of view, they coincide exactly with the spin— (j = ^) representation of 1(70(1 together with 

c = 1 , T b — ► -T . 

This should not be surprising as the symmetry between the holomorphic/antiholomorphic sector must be restored 
at the boundary. Due to the presence of the boundary which breaks the U qo (sl2) symmetry, the quantum group 
symmetry of the boundary sine-Gordon model with fixed boundary conditions is generated by the non-local 
conserved charges (|70|) for the special case (symmetric) cq — 1. In particular, single representations are associated 
with a restricted Askey- Wilson quadratic algebra. As before, the non-local boundary conserved charges satisfy 
the q— deformed Dolan-Grady relations JHJ) with ((69(1 and cq = 1. Conservation of the quantities G^ 2r ' for cq = 1 
implies that the boundary space representation has dimension 2j + 1. A convenient basis for the boundary 
(degenerate) ground state |0)s is provided by the Askey- Wilson polynomials with discrete weights, or can be 
written in the same manner than 1)75(1 . 

The boundary scattering properties are known to be encoded in the boundary reflection matrix, which can be 
obtained using the quantum reflection symmetry |3f)j (see also |32j). In the formalism of Section 3, it corresponds 
to the choice of representation tt^ 1 ^ 2 ^ x in 144(1 where tt^(A) ~ e+, n^°\A*) ~ e_ for the non-local charges. 
Using the convenient parameterization of [53j , it is straightforward to relate ((13(1 with the reflection matrix of |16| 
using the appropriate unitarization factor. 

• Massless limit: anisotropic Kondo and Bazhanov-Lukyanov-Zamolodchikov models. The mass- 
less limit of the model ((78(1 can be reached in various ways, according to the realization of the boundary opera- 
tors 18 one would like to have. In the dynamical boundary case, it should be reminded that the integrability is 
preserved for boundary operators which satisfy asymptotically ((20(1 and the substitution A — > £ + , A* — > £_ with 
p = p* = c?{q x / 2 — q -1 / 2 ) (see 5B) for details). The "massless" limit of ((78(1 is obtained if the eigenvalues of 
fi 1 ^ 2 £±\ l ji-,o are finite. This can be done, for instance, by taking the limit vq — > oo in ((68(1 such that VqP ~ 1. It 
gives 

A c+v Q S + q S3/2 and A* -> c_v S*_q" S3/2 (82) 
where ((45(1 . Using ((82(1 in ((78(1 . the resulting Hamiltonian becomes 

H massless ~ -L f dx((n(x)) 2 + (d x ^x)) 2 ) ^ 2 {S + q^ 2 e^/ 2 + S _ g -«M^) (83 ) 

For 9=1, this Hamiltonian is the bozonized version of the anisotropic Kondo model |54| . For q ^ 1, one obtains 
the spin— j generalization studied in details in (55j . 

An other example of massless model is obtained using a different realization of the Lax operator in ll'-SI) , and 
proceeding by analogy to obtain 1(65(1 . Taking the appropriate limit, one obtains 

A - i{q 1 / 2 - g-VajVV and A* -> i(q^ 2 - q-^^a.q^ 2 , (84) 

where the operators a + , a_ and q N generate the q— oscillator algebra 

a_a + — <7 ±:L a + a_ = q^ N and q N a±q~ N = q ±1 a± ■ (85) 

Notice that q x / 2 AA* — q~ x l 2 A* A = q 1 ! 2 — q~ x / 2 . For this realization of the boundary operators (in which case (0 
simplify to q— Serre relations), the model 1(78(1 reduces to the one proposed by Bazhanov-Lukyanov-Zamolodchikov 
in (2()j at zero voltage 19 . This model finds interesting application in non-equilibrium systems. 

17 These parameters are expressed in terms of the boundary parameters in the Lagrangian and some bulk contributions. I thank Z. 
Bajnok who attracted my attention to this point. 

18 I thank H. Saleur for stimulating communications about this problem. 

19 In case of non-zero voltage, one can use a slightly different expression for the generating function. 
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• Reflectionless points: Bassi-LeClair massive Kondo model. The reflectionless points correspond to 
the values of the coupling constant (3 is such that q = ±1. In |56| a rather general massive version of the anisotropic 
spin 1/2 Kondo model has been proposed. At the reflectionless points, it was claimed to be integrable unlike the 
massless version which remains integrable for arbitrary values of f3. Using (|65J) . one can see that the Hamiltonian 
associated with (178(1 coincides with the model proposed by Bassi and LeClair in j^H] for q = 1. A more interesting 
property can furthermore be exhibited: at this value of q, the boundary operators £± and boundary non-local 
conserved charges Qi"* generate the infinite dimensional Onsager algebra (0- It is then possible to construct the 
corresponding conserved quantities in involution which ensure the integrability of the theory and confirms the 
conjecture of (SSj- Furthermore, this might open a way to extract exact information (correlation functions,...) in 
the massive regime. 

• Remark on the decoupling limit: Bulk/Boundary "dual" description. It is rather instructive to 
consider the limit in which the bulk degrees of freedom (sine-Gordon field) and the boundary ones £±(y) decouple. 
This can be realized either for Neumann boundary conditions, or vanishing coupling constant (3 — in (|78|l . For 
Neumann boundary conditions, i.e. £±(y) = 0, only the bulk sine-Gordon contribution and the term Aboundary 
survive. The non-local conserved charges in this limit (35] simplify to Q±} \g ± _ tConst in (|81|) . On the other hand, 
in the second case j3 = 0, the action 1(78(1 becomes 

/oo 
dy(£-(y) + £ + (y)) + Aboundary ■ (86) 
-oo 

In the asymptotic limit {x,y) — > ±00, it is easy to check that £± (for j3 = 0) and Q± \s+—^Const (for Neumann 
boundary conditions) satisfy exactly the same algebraic relations Q. In view of the weak-strong coupling duality, 
in the perturbative regime /3 — > in 178(1 , one can see that the boundary contribution 1(86(1 actually contains some 
informations about the nonperturbative regime which is usually characterized by the non-local conserved charges. 
This might explain why, for the analytic continuation (3 — > ib in l|78|l . the dynamical boundary sinh-Gordon model 
is self-dual under weak-strong coupling duality (the reflection amplitude of the fundamental particle is invariant 
under the change b <-> 2/6) ^H], contrary to the model with fixed boundary conditions \b'6\ . 



6 Concluding remarks 

The spectral parameter dependent reflection equation is usually considered in the context of integrable systems 
with boundaries. However, based on the analysis of this quadratic algebra we have exhibited a new quantum 
integrable structure which is not specific to boundary systems. The integrability condition was shown to be 
associated with q— deformed Dolan-Grady relations J5j. Our construction, which extends the one proposed by 
Dolan-Grady some years ago [5], leads to the existence of an (in)finite set of conserved quantities in involution 
40 .In particular, it shows that the concept of "superintegrability" can be further extended for q ^ 1. In this 
direction, it is a challenge to find the corresponding (in)finite dimensional q— Onsager algebra that generate these 
conserved quantities in massive theories, and would provide an alternative approach to massive quantum integrable 
models. This new hidden symmetry JSJ) was exhibited in various bulk or boundary quantum integrable models 
in which case the operators A, A* have been identified and the first conserved quantities l@J and 1(52(1 constructed 
explicitly. Among the interesting applications, we have shown that the sine-Gordon model enjoys a remarkable 
dynamical symmetry associated with the quadratic Askey- Wilson algebra. From this, we have obtained the 
asymptotic one-particle states in terms of q— orthogonal polynomials, as well as the structure of the boundary 
space of state in this model with a boundary. From a more general point of view, we have also shown that several 
examples of massive/massless boundary integrable models are special cases of the model 1(78(1 . Finally, we have 
pointed out a bulk-boundary "dual" relation in the sine-Gordon model with a dynamical boundary. Although we 
didn't discuss the XXZ spin chain with/without a boundary here, a similar description also holds in this model. 

An interesting problem is wether the present approach can be generalized to higher rank affine Lie algebras. 
Actually, there has been some interest in generalizations of the Onsager algebra and corresponding Dolan-Grady 
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relations. Indeed, the Onsager algebra is intimately related with 5/2- In case of sl ni such extension was studied 
in details in |58| . There, an integrable Hamiltonian of the form 

H = K4A1 + k 2 A 2 + ... + k„A„ , n > 3 (87) 

was proposed where Ki are arbitrary constants. This Hamiltonian is a member of an infinite family of commuting 
integrals of motion if the operators A^ satisfy certain trilinear relations. From our results, it is tempting to consider 
a higher rank q— deformed generalization of JSJ with the following trilinear q— deformed relations: 

[AiJAi.Aj],],-! = piAj , and [Ai,Aj]=0 if |i-j|>l (88) 

where {i,j} £ {1, ...,n}. In some sense, these relations extend the concept of Leonard pair to any finite Lie 
algebra. Obviously, using an appropriate shift A, — > A, + a, it is not difficult to generalize the relations 1)20(1 . 
For the same reasons as before, it is consequently sufficient to focus on the relations 1(88(1 . For q = 1 in 1(88(1 . 
one obtains the relations found in J^Hj- Also, for pi — 0, one recovers the q— Serre relations associated with 
U q (sl n ). We intend to consider this model together with applications to massive Toda field theories in a separate 
publication. 

To conclude, we would like to stress that the power of the relations (JSJ) lie in the algebraic statement which 
does not refer to the number of dimensions or the the nature of the space-time (continuous, discrete,...) "hidden" 
behind the operators A, A*. An interesting problem is consequently to find some examples where such hidden 
symmetry appears explicitly in higher dimensions. 

Note added: 

Using the realization l(54|) . it is easy to relate our general expression ((13(1 with ((22(1 to known results. For instance, 
choosing Co = TL c i = T2/(<7 — q^ 1 ), c 2 = and q — > q~ 2 it is not difficult to check that lf^)l coincides with the 
one proposed in [55] (for (— )) (see also Q7]) or the one that follows from I84 (for (+)). 

Also, it should also be noted that the reflection matrix 1(13(1 . for some special realizations of Leonard pairs in 
terms of Uqi/z^sfa) can be related with a solution proposed some years ago by A. Zabrodin in 49 . This solution 
was obtained from the non-dynamical K— matrix in |16| and L-operators satisfying ((32(1 . It should be stressed 
that our solution ((13(1 is derived independently of a realization of Leonard pairs, either directly from the RE ((12(1 
or using the quantum affine reflection symmetry described in Section 3. 

Acknowledgements: I am grateful to V. Bazhanov for discussions at the preliminary stage of this work and 
P. Terwilliger for explanations about his work on tridiagonal algebras and Leonard pairs, as well as his advice 
in improving the first version. I also thank S. Nicolis for explanations about the Azbel-Hofstadter model, M. 
Niedermaier, P. Dorey and especially P. Forgacs for discussions, interest in this work and important comments 
about the manuscript. Part of this work is supported by the TMR Network EUCLID "Integrable models and 
applications: from strings to condensed matter", contract number HPRN-CT-2002-00325. 

Appendix 

In this appendix, we want to show that an Hamitonian of the form ||3J can be obtained starting from an 
L— operator instead of a if— operator. Indeed, in the quantum inverse scattering method, an integrable lattice 
theory can be constructed considering the monodromy matrix 

T{u) = L VoN (u)...L Vo2 (u)L Vq1 (u) (89) 

where L Vg j (it) acting on the quantum space Vj is a solution of the Yang-Baxter algebra JHJ for Vo = Vq and 
V — > Vj. Taking the trace over the auxiliary space Vo of 1(89(1 . one obtains the transfer matrix 

T(u)=Tr Va (T(u)) (90) 
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acting on ®^ = iVj. This generating function provides a set of mutually commuting conserved quantities. For 
instance, let us focus on the Yang-Baxter algebra © with (J7J) associated with [^1/2(5/2) in the spin— ^ auxiliary 
space representation Vo = Vq and N = 1 in (|89|l . Following 03], one can consider 

W \ u fcg 1 /4( P g-i)-i/2_ u -i fc -i (7 -i/4 ( pQ-i ) i/2 -inp l l 2 Q- 1 -i^p^p- 1 ) 1 ' 

with parameters {£;, k, k*} 6 (F. Using the realization of the Leonard pair A, A* in terms of the Weyl algebra W q 
l|54|) , one obtains from (|90|) a Hamiltonian acting on V of the form 

H=(q-q- 1 )( K A + K*A*) (92) 

where A, A* satisfy the q— deformed Dolan-Grady relations 0. 
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